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ABSTRACT 



^> In this paper we present a model of the /3 Leo debris disc, with an emphasis on modelUng 

r^ the resolved PACS images obtained as part of the Herschel key programme DEBRIS. We 

also present new SPIRE images of the disc at 250 \ini, as well as new constraints on the disc 

1-^ from SCUBA-2, mid-IR and scattered light imaging. Combining all available observational 

*~-i constraints, we find three possible models for the (3 Leo (HD 102647) debris disc: (i) A 2 com- 

Q ponent model, comprised of a hot component at 2 AU and a cold component from 15-70 AU. 

?-H (ii) A 3 component model with hot dust at 2 AU, warm dust at 9 AU, and a cold component 

c/3 from 30-70 AU, is equally valid since the cold emission is not resolved within 30 AU. (iii) A 

Ci somewhat less likely possibility is that the system consists of a single very eccentric planetes- 
imal population, with pericentres at 2 AU and apocentres at 65 AU. Thus, despite the wealth 

p^ of observational constraints significant ambiguities remain; deep mid-IR and scattered light 

^ imaging of the dust distribution within 30 AU seems the most promising method to resolve 

\j^ the degeneracy. We discuss the implications for the possible planetary system architecture; 

T— I e.g., the 2 component model suggests planets may exist at 2-15 AU, while the 3 component 

C^ model suggests planets between 2-30 AU with a stable region containing the dust belt at 9 

^D AU, and there should be no planets between 2-65 AU in the eccentric planetesimal model. 

r — We suggest that the hot dust may originate in the disintegration of comets scattered in from 

f^ the cold disc, and examine all A stars known to harbour both hot and cold dust to consider 

7—i the possibility that the ratio of hot and cold dust luminosities is indicative of the intervening 

^-H planetary system architecture. 

^ 1 INTRODUCTION of structure and evolution of planetary systems. Resolved images 



X 



can indicate that infrared excess is being produced by multiple dust 
Debris discs are distributions of dust and planetesimals with radii populations and can also break the degeneracy between the radial 
^ of 1-1000 AU around main sequence stars (See Wyatt 2008 for a i^^^ji^jj ^^ ^^^ ^^^^ ^^^ ■^^^ temperature, 
recent review). The dust grains in these discs are small, and so can- 
not be primordial as they would have been blown out of the system The DEBRIS (Disc Emission via a Bias-free Reconnaissance 
by radiation pressure on timescales shorter than the stellar age. The ^ *e Infrared/Sub-mm) survey (Matthews et al. in Prep, Phillips 
dust must be continually replenished from a population of colliding et al. 2010), is an Open Time Key Program on the Herschel Space 
planetesimals, thought to contain bodies up to ~ 1km in size (Wyatt Observatory which uses PACS (Photodetector Array Camera and 
& Dent 2002). However, as debris discs age the planetesimal pop- Spectrometer Poglitsch et al. 2010) and SPIRE (Spectral and Pho- 
ulation is ground down, so discs become less massive and fainter tometric Imaging REceiver, Griffin et al. 2010) to detect, resolve 
(Dominik & Decin 2003). At far-infrared and sub-millimetre wave- and characterise debris discs around a volume-limited sample of 
lengths debris discs are optically thin, the disc to star contrast is 446 A through M stars (Matthews et al. 2010). /3 Leo (HD 102647) 
favourable, and these wavelengths are sensitive to the large (up to ^^s observed at 100 ^m and 160 tim with PACS as part of the DE- 
~ 1mm) grains that dominate the dust mass in debris discs. BRIS survey Science Demonstration Phase (Matthews et al. 2010). 
The dust morphology of a debris disc can be shaped by planets /? Leo (A3V, L«=14.OL0) is a 5 Scuti type star at a distance 
in the system so resolved images of discs help constrain models of 11.1 pc. The infrared excess around this main-sequence star 
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was first discovered by IRAS (Infrared Astronomical Satellite Oud- 
maijer et al. 1992), then confirmed with Spitzer (Su et al. 2006). 
The excess was unresolved in mid-IR imaging (Jayawardhana et al. 
2001; Akeson et al. 2009, § 2.2) although differences between the 
IRAS and ISO (Infrared Space Observatory) fluxes led Laureijs 
et al. (2002) to suggest that the disc emission may be somewhat 
extended in the ISO beam (52" aperture). Chen et al. (2006) ob- 
tained a Spitzer IRS (Spitzer Infrared Spectrograph) spectrum of /? 
Leo (See § 3.2) and found a featureless continuum spectrum con- 
sistent with dust at ~120 K located at 19 AU from the star. A very 
hot excess has also been partially resolved using infrared interfer- 
ometry with the FLUOR (Fiber Linked Unit for Recombination) 
instrument at the CHARA (Center for High Angular Resolution 
Astronomy) array at 2 (i.m and BLINC (Bracewell Infrared Nulling 
Cryostat) at 10 ^m (Akeson et al. 2009, Stock et al. (2010), see § 
2.8 for more details). 

P Leo is thought to be a member of the IC2391 moving group, 
giving an age of 45 Myr (Nakajima et al. 2010). The age of this 
source determined from isochrone fitting is 50-331 Myr (Lachaume 
et al. 1999; Song et al. 2001). Di Folco et al. (2004) derive an age 
from measuring the stellar radius and suggest 100 Myr. We as- 
sume an age of 45 Myr for this paper. (3 Leo does not have any 
known companions within a few arcsecs of the star. The Wash- 
ington Double Star (WDS) catalogue lists three companions with 
common proper motion for /3 Leo. These stars are located from 
40" to 240" from the primary with V magnitude differences of 
6.3 to 13 (Worley & Douglass 1997); these stars are not, however, 
physically associated with /? Leo (Phillips et al. 2010). /? Leo was 
also included in a high precision radial velocity survey of early type 
dwarfs with HARPS (Lagrange et al. 2009) for planetary or brown 
dwarf mass candidates and was found to have no companions with 
mass >4.2Mj„p with periods <10 days with 99.7 percent probabil- 
ity. 

This paper uses multi-wavelength modelling of the /3 Leo de- 
bris disc including PACS and SPIRE observations to present a self- 
consistent explanation of the system. In §2 we present the observa- 
tions, including 100 and 160 \im PACS imaging (previously pub- 
lished in Matthews et al. 2010) and archive Gemini MICHELLE 12 
\im and 18 \im imaging. In §3 we confront the observations with 
models to determine the disc parameters. In §4 we discuss the im- 
plication of the inferred structure for the status of planet formation 
in this system. 



2 OBSERVATIONS 

2.1 Fitting the Stellar Photosphere 

Being the 5th nearest A-type star, /3 Leo is bright and saturated in 
modern surveys such as 2MASS. In addition, the infrared excess 
has a hot component, which contributes to bands typically used for 
modelling stellar photospheres (Akeson et al. 2009; Stock et al. 
2010). We use the equivalent 2MASS Ks magnitude of 1.93 from 
Stock et al. (2010), and include mean UBV and Hipparcos Hp pho- 
tometry (Mermilliod 2006; ESA 1997). The best fitting Castelli & 
Kurucz (2003) model, found by a x^ minimisation method, has 
Teff = 8660K, L-, = 14Lq, and i?* = 1.66i?Q. In the MIPS 24 
^m band we predict 1171 ± 15 mjy, 1% lower, but consistent with 
the value in Stock et al. (2010). The predicted stellar flux densities 
at the PACS effective wavelengths of 100 and 160 /im are 64 ± 1 
and 26 ± 0.5 mJy respectively. These are 9% higher than the val- 
ues presented in Matthews et al. (2010) due to improvements in the 
fitting procedure and data used. 



2.2 Herschel PACS Observations 

Observations of the /3 Leo disc at 100 (i.m and 160 |J.m were taken 
with the ESA Herschel Space Observatory (Pilbratt et al. 2010) 
using the PACS (Poglitsch et al. 2010) instrument in photometry 
mode as part of the Science Demonstration Phase observations for 
the DEBRIS survey. DEBRIS is a flux-limited survey which ob- 
serves each target to a uniform depth of 1.5 mJy/beam at 100 |J.m. 
The images were first presented in Matthews et al. 2010. /? Leo 
was observed in small scan-map mode (see PACS Observers' Man- 
ual^ for details). Scan map observations had eight repeats in a sin- 
gle scan direction at a scan rate of 20"/s. Four 3' scan legs were 
performed per map with a 2" separation between legs. The total 
observing time was 1220s. 

These data were reduced using the Herschel Interactive Pro- 
cessing Environment (HIPE Ott (2010)). Maps were obtained via 
the default PACS naive map-making method photProject in HIPE. 
The data were pre-filtered to remove low-frequency (1/f) noise us- 
ing a boxcar filter with a width of 98". All bright sources in the map 
were masked prior to filtering to avoid ringing type artifacts. 

Photometry on the maps of /3 Leo (Figure 1) using a 13" ra- 
dius circular aperture centred on the emission peaks yielded fluxes 
of 480±30 mJy at 100 Hm and 215±32 mJy at 160 ^m, respec- 
tively. These fluxes greatly exceed the rms noise levels in the maps 
of 1.4 miy/beam at 100 (i.m and 3.1 mJy/beam at 160 |J.m. These 
values have been colour corrected using the values from the PACS 
Observers Manual assuming a temperature of 120K corresponding 
to the temperature of the disc estimated from a blackbody fit to 
the SED. The quoted errors on these fluxes do not include calibra- 
tion uncertainties which are estimated to be 10 and 20% at 100 and 
160 (im (Poglitsch et al. 2010). These were combined in quadrature 
with statistical uncertainties from the rms levels in the maps. The 
stellar flux for fitting the photosphere is 64 mJy and 26 mJy at 100 
and 160 (i.m, respectively. This gives excess fluxes of 416±30 mJy 
at 100 (im and I89±32 mJy at 160 (J-m. These values differ from 
those in Matthews et al. 2010, which quoted 500±50 mJy at 100 
(im and 230±48 mJy at 160 (im in a 20"radius aperture. This dif- 
ference is because the data presented in this paper were rereduced 
using a different filter scale and a more recent version of HIPE 
(4.2.0 developer's build) and a smaller aperture was used. 

The maps of /? Leo (Figure 1) appear extended compared to 
the beam (indicated by the circle in the top left corner of the im- 
ages). At 100 M-m the nominal beam size is 6. "6 x 6. "9. At 160 
(j.m it is 10. "7 x 12." 1. We have obtained 2 bright point source 
images (Vesta and a Boo) taken in the same observing mode as the 
P Leo data to serve as PSF references. These have beams sizes of 
6. "6 X 6. "9 at 100 [im and 10. "5 x 11.5" at 160 pun for Vesta and 
6. "5 X 6. "8 at 100 t^m and 10. "2 x 11. "7 at 160 ^m for a Boo. 
Fitting a Gaussian to the /? Leo image gives a FWHM at 100 |im of 
9.2"±0.1 X 10.4"±0.1. At 160 pun the FWHM of a Gaussian fitted 
to the image is 13.6"±0.2x 12.0"±0.2. Figure I show the /? Leo 
images after subtraction of the a Boo PSF scaled to the appropriate 
stellar flux. The PSF has an asymmetric three lobed structure and 
was rotated to the same spacecraft angle at which the /? Leo data 
were taken at to ensure accurate subtraction. 

/3 Leo appears extended in all directions with respect to the 
PSF, with the major axis of the disc at a position angle of 125° ±15° 
E of N, which is the mean of the position angles of the Gaussians 
fitted at 100 (PA: 118°) and 160 ^m (PA: 132.5°). The ratios of the 

-'- PACS operating Manual:http : //herschel . esac .esa. int/ 

Docs /PACS /html /pacs_om. html 
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Figure 1. Images of the 100 (left) and 160 |J.m (right) emission for /3 Leo taken with PACS. The images have had a PSF scaled to the stellar flux (64 mJy and 
26 mJy at 100 and 160 |J.m) subtracted, and hence are maps of the excess emission. The pixel scale in 1" per pixel at 100 |J.m and 2"per pixel at 160 |j.m. 
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semi-minor to semi-major axis of the Gaussians (0.84 at 100 i^m 
and 0.88 at 160iJ.m) give a mean inclination of 57±7° from edge 
on. The stellar subtraction was also done using the Vesta PSF to 
examine the effect of PSF variation on the residuals, and resulted 
in no significant change in the width of a Gaussian fitted to the disc, 
and no significant change in the flux in a 13"radius aperture. 



2.3 MICHELLE Mid-Infrared Observations 

Mid-IR observations of 13 Leo taken with MICHELLE on Gem- 
ini North for /? Leo were retrieved from the Gemini Science 
Archive (Program Names: GN-2007A-C-10 [PI: Beichman] and 
GN-2006A-Q-10 [PI: Moerchen]). These data are detailed in Ta- 
ble 1 and were taken with filters Qa (Ac=18.1 vim, AA=1.51 |i.m) 
andN' (Ac=lL3 Hm, AA=1.07 ^.m). MICHELLE has a pixel scale 
of 0.1005" per pixel. Imaging was taken using a ABBA chop-nod 
sequence with a position angle of 30° for all observations. The 
data taken under GN-2006A-Q-10 were divided into 3 observa- 
tion groups, with observations of a standard star from Cohen et al. 
(1999), HD 98118 (MOIII, Fi8^m=4.2 Jy), taken before and after 
the observations of /? Leo to serve as a standard star for flux calibra- 
tion and to monitor the PSF. The data taken under GN-2007A-C- 10 
were previously published by Akeson et al. (2009). Observations of 
HD109511 (K0,Fi8^m = 1.4Jy) preceded and followed the obser- 
vations of /3 Leo as a standard star and PSF reference. Combining 
all observations, the total on source time at 18 vim is 1751s and 
1175s at 12 ^m. 

The data were reduced using custom routines described in 
Smith et al. (2008). The data reduction involved determining a gain 
map for each observation using the mean values of each frame to 
construct a map of pixel responsivity. The on-source pixels were 
masked during this process, making this equivalent to a sky flat 
field frame. A DC offset was then determined by calculating the 
mean pixel values in every row and every column, again mask- 
ing pixels where there was source emission present. This was then 
subtracted from the final image to ensure a flat background. Pix- 



els which showed high or low gain in comparison with the median 
response throughout the observation were masked off. To avoid er- 
rors in co-adding the data which could arise from misalignment of 
the images, we fitted a Gaussian with a sub-pixel centroid to accu- 
rately determine the centre of the image and so the position of the 
star. The re-binning was done using bilinear interpolation across 
the array. 

The observations were divided into four groups (as shown in 
Table 1) such that flux calibration for each group was done with 
standards observed at a similar airmass to P Leo. The calibration 
levels were compared with the airmass for the standard star obser- 
vations and no correlation was found, so no extinction correction 
was applied to the calibration factors. 

For the images of /3 Leo in each group a calibration factor was 
determined using a co-add of the two standard star observations in 
that group. The average calibration error for the 3 18 vim groups 
was 8 %, but as seen in Table 1 the large calibration uncertainties 
lead to a wide range in fluxes among the three groups. The 12 vim 
group had a calibration error of 1 1%. These centred, flux calibrated 
images from each group were then used to produce final images of 
P Leo and the standard star. The first standard (Stdl) and image 
(Iml) in Group 1 (See Table 1) at 18 ^m and Stdl2.1 and Iml2.1 
at 12 vim showed elongation in the telescope chop-nod direction, 
so were not used in the final co-added images, but were used when 
calculating the flux calibration factors. 

Photometry was performed on the final co-added images of 
P Leo using a 1.0 arcsec radius circular aperture centred on the 
star. The stellar flux is expected to be 5297 mJy and 2020 mJy at 
1 1.3 and 18.1 vim respectively, from fitting a model spectrum as de- 
scribed in § 4. Our photometry yields fluxes of 5822±476 mJy at 
at 12 vim and 2360±3 12 mJy at 18.6 vim including both calibration 
and photometric errors, giving an excess of 525±476 mJy at 1 2 vim 
and 340±312 mJy at 18 vim- The statistical noise was determined 
using an annulus with an inner radius of 7" and an outer radius of 
8" centred on the star, resulting in an error averaged over all the 
groups of 0.31 mjy/arcsec^at 12 vim and 0.41 mjy/arcsec^ at 18 
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Table 1. Observations taken under proposal GN-2006A-Q-10 and GN-2007A-C-10 in order. The integration time listed is the on-source integration time. 
Fluxes are for a 1" radius aperture centred on the star The group indicates the standard-science-standard observing pattern used. 



Program 


Date 


Object 


Group 


Name 


Filter 


Integration Time (s) 


Calibrated Flux (Jy) 


GN-2006A-Q-10 


10/05/06 


HD 98118 


1 


Stdl 


QA 


82 


4.20 


GN-2006A-Q-10 


10/05/06 


/3Leo 


1 


Iml 


QA 


396 


3.17 ±0.13 


GN-2006A-Q-10 


10/05/06 


/3Leo 


1 


Im2 


QA 


245 


3.17 ±0.13 


GN-2006A-Q-10 


10/05/06 


HD98118 


1 


Std2 


QA 


82 


4.20 


GN-2006A-Q-10 


15/05/06 


HD98118 


2 


Std3 


QA 


82 


4.20 


GN-2006A-Q-10 


15/05/06 


/3Leo 


2 


Im3 


QA 


326 


2.80±0.25 


GN-2006A-Q-10 


15/05/06 


HD98118 


2 


Std4 


QA 


82 


4.20 


GN-2006A-Q-10 


11/06/06 


HD98118 


3 


Stdl2.1 


N' 


47 


3.65 


GN-2006A-Q-10 


11/06/06 


/3Leo 


3 


Iml2.1 


N' 


376 


5.82±0.51 


GN-2006A-Q-10 


1 1/06/06 


/3Leo 


3 


lml2.2 


N' 


376 


5.82±0.51 


GN-2006A-Q-10 


1 1/06/06 


/3Leo 


3 


lml2.3 


N' 


423 


5.82±0.51 


GN-2006A-Q-10 


11/06/06 


HD98118 


3 


Stdl2.2 


N' 


47 


3.65 


GN-2007A-C-10 


8/03/07 


HD109511 


4 


Std5 


QA 


86 


1.4 


GN-2007A-C-10 


8/03/07 


/3Leo 


4 


Im4 


QA 


392 


2.36±0.33 


GN-2007A-C-10 


8/03/07 


/3Leo 


4 


Im5 


QA 


392 


2.36±0.33 


GN-2007A-C-10 


8/03/07 


HD98118 


4 


Std6 


QA 


86 


1.4 



Hm. The IRS spectrum presented in Stock et al. (2010) gives an ex- 
cess above the photosphere of 61±103 mjy at 12 \un and 256±53 
mjy at 18 |J.m. The photospheric fit used here is consistent with that 
used in Stock et al. (2010). The MICHELLE photometry presented 
here therefore agrees with the IRS results, but the calibration errors 
are too large to detect the excess. To assess whether these mid-IR 
images have resolved the disc, line cuts were taken along the PA 
of the extension (125°E of N) seen in the PACS images (125°) for 
both /3 Leo and the standard star . These are shown in Figure 2. fi 
Leo shows no extension when compared with the PSF, and the ex- 
cess at both 12 and 18 |J.m is unresolved. This is consistent with the 
previous results of Akeson et al. (2009); Moerchen et al. (2010). 
None of the groups show significant extension. We also tried sub- 
tracting the standard star image scaled to the peak flux of the /3 Leo 
images, but no significant structure remains. 

As the location of the disc is known from the resolved PACS 
images with peak emission at ~5" radius, we tried convolving the 
co-added images with a series of Gaussians with FWHMs from 1" 
to 20" to find any large scale, low surface brightness features in the 
outer regions of the MICHELLE images. The MICHELLE field 
of view is 32"x24", but the chop throw of 15" limits the usable 
region. We find no coherent features at the expected location of the 
debris disc in the convolved image. The surface brightness in an 
annulus from 0.5" to I" in the MICHELLE 18 (i.m imaging of P 
Leo is 0.41 ±0.34 mjy/pixel. However, the brightness in the same 
annulus of the PSF scaled to the flux of /3 Leo is 0.3 1 mjy/pixel. 



2.4 SCUBA 2 

The SCUBA 2 instrument (Sub-mm Common User Bolometer Ar- 
ray 2) (HoUand et al. 2006) operating on the JCMT ((James Clerk 
Maxwell Telescope) was used to obtain fully sampled maps of /? 
Leo, and the surrounding region at 450 /.im and 850 /im. These data 
were obtained on the nights of the 1 6th and 1 7th February 2010 as 
part of the shared risk observing phase of instrument commission- 
ing. The final maps contain data from six observations obtained on 
each night, giving a median integration time of 3 hours per pixel 
in the central one square arcminute, and a total observing time of 6 
hours. The 'daisy' pattern observing mode was used at a scanning 
rate of 120 arcseconds per second (Kackley et al. 2010), giving a 
map with a usable area of ~9 square arcminutes. 



Figure 2. The profile of the line cuts through the total co-added Mid-IR im- 
ages at 125°, the position angle of the extension seen in the PACS images. 
The dotted lines are the final co-added image of /3 Leo in each group and 
the solid lines are the final coadded standard star in each group. The 1 2 |j.m 
line cut is the top image, the 18 |J.m is the bottom image. /3 Leo shows no 
significant extension with respect to the PSF, indicating the excess at this 
wavelength is unresolved. In fact, the observations appear narrower than 
the PSF, but this is not significant and is due to variation in observing con- 
ditions. 
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The data were processed and calibrated using the SMURF 
package in the Namaka Starlink release (Jenness et al. (2010); 
Dempsey et al. (2010)). The data were high-pass filtered to miti- 
gate the 1/f noise present in the data, with filter parameters set so 
as to retain features smaller than ~I20 arcseconds. Maps were ob- 
tained at 450 [im and 850 (xm and the data were reduced using the 
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Figure 3. Spitzer MIPS observations of /3 Leo at 24 [iva. The disc is detected and partially resolved compai'ed to the PSF. The disc image is shown in the left 
panel, the PSF is the centre panel and the lower panel shows the surface brightness profiles of the Observations (solid line) and the PSF (dashed line) showing 
the filling in of the first dark Airy ring between 5.5"and 7.5 ". 




RA Offset, arcsec 



RA Offset, arcsec 



-2.2e+01 (linear} 1.1e+04 

(a) Spitzer MIPS 24 |xm Observations of /3 Leo. 
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(c) Surface brightness profile of 24i^m /3 Leo observation 
and the PSF 



latest version of the SMURF pipeline. /3 Leo was not detected in 
either band. Photometry on the 450 \im map gives a 3 ct upper limit 
of 50 mjy/beam At 850 \un the 3 a upper limit is 6 mJy/beam. 
The expected stellar flux in these wavebands is 3.1 mjy and 0.8 
mjy respectively. Holmes et al. (2003) found an upper limit of 20 
mJy at 870|i.m, so these upper limits significantly improve the con- 
straints on the SED in the sub-mm and were used in the SED fitting 
described in §3. 



2.5 Spitzer MIPS Observations 

In addition to the photometric points at 24 and 70 |i.m on the SED 
shown in Figure 1 , we also include the partially resolved MIPS im- 
age of the disc at 24 |i.m, presented in Stock et al. (2010). Stock 
et al. (2010) observed that the 24 \un photometric point obtained 
using an aperture of 14". 94 gives a total integrated flux of 1623±33 



mJy which is 2.5% higher than the flux obtained in a 6."23 aper- 
ture, leading to the suggestion that the disc may be marginally re- 
solved at this wavelength. The photometric value from the larger 
aperture is used for SED fitting. The MIPS 24 \un image is shown 
in the first panel of Figure 3. Compared to the PSF shown in the 
right panel of Figure 3 the first dark Airy ring (between radii of 
5.5"to 7.5") appears more filled in. Subtraction of the PSF scaled 
to the expected photospheric flux of /3 Leo gives residuals with a 
FWHM of 6.88" x6.6l"at a position angle of 1 18° which is larger 
than the FWHM of the PSF (5.69" x5. 53"). This is again consis- 
tent with a disc that is marginally extended with an inclination of 
60±10° from edge on and a position angle of 1 18±3°. Comparing 
P Leo to the PSF observed for a star with hot dust at <2 AU (i.e. 
similar to the inteferometric detection) of the same spectral type ((^ 
Lep - 5."61 x5."55 Su et al. 2008) and to a theoretical PSF (gen- 
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erated using Tiny Tim ^) also suggests that the disc is marginally 
resolved. If the 24 |j.m observations are deconvolved from the PSF 
(i.e. the FWHM of the PSF is subtracted in quadrature from the 
FWHM of the observations), this gives the FWHM of the residuals 
as 3.9". The position angle and inclination of these residuals are 
consistent with the disc resolved using PACS. The 70 [xm Spitzer 
MIPS obsevations are consistent with a point source. We include 
the 24 fi.m image when considering the observational constraints 
for a model of the B Leo debris disc. 



2.6 Optical Observations 

/3 Leo was imaged using the Hubble Space Telescope Advanced 
Camera for Surveys High Resolution Channel (HRC) coronagraph 
on March 25, 2004 (GO-9475, PI Kalas). We occulted /3 Leo with 
the 1.8" diameter occulting spot approximately centered on the 
1024x1024 pixel camera. Seven 140 second exposures were ac- 
quired using the F606W filter (Ac = 591 nm, <5A = 234 nm). 
The AlV star HD95418 was observed before j3 Leo to serve as a 
reference for the stellar PSF . Cosmic rays were filtered by taking 
the median combination of the seven /3 Leo exposures, as well as 
the six 180 second exposures on HD95418. It should be noted that 
HD95418 /? UMa has an IR-excess indicative of a debris disc but 
this is not detected in scattered light. The HD 95418 observation 
was registered and scaled by a factor of 1.22 to subtract the /3 Leo 
PSF. The resulting, PSF-subtracted image was then corrected for 
geometric distortion giving a pixel scale of 0.025" pix^^. Figure 5 
shows the resulting optical image of /3 Leo. We find a residual halo 
of light in an annulus 2.5 — 4.0" from the star. This halo could plau- 
sibly originate from either dust scattered light or an imperfect PSF 
subtraction. The latter effect results from thermal "breathing" of 
the telescope between the observations of the two stars, as well as a 
small color mismatch between /3 Leo and HD 95418. The limits on 
surface brightness from these observations were used as constraints 
for the Modelling described in §3. 



2.7 Herschel SPIRE Observations 

Photometric observations of j3 Leo were obtained with the SPIRE 
instrument on Herschel (Griffin et al. 2010), providing maps in 
wavebands centred at 250, 350, and 500 /xm. These data were 
observed on the 23rd November 2011 (Herschel operational day 
558) using the 'small map' observing mode. Five repeat maps were 
performed, effectively giving confusion noise limited maps in all 
bands. The data were reduced using HIPE Version 6.0, build num- 
ber 1985, and the standard pipeline script was used to perform 
the reduction. The images are shown in Figure 5. Photometry was 
performed on the images. At 250 \im a 14"radius circular aper- 
ture (beam size 18.9"x 17.6"^) gives 51±12 mjy. At 350 ^m a 
30" aperture (beam size 25.6" x 24.2") gives I8±7 mJy which gives 
a 3(T upper limit of 39 miy. At 500 (xm a 40"aperture (beam size 
38.0" X 34.6") gives 2.8±4. 1 mJy, which corresponds to a 3a upper 
limit of 15.1 mJy. 

To assess whether the 250 (im image had resolved the disc 
we fit the image with a 2D Gaussian, which has a FWHM of 



^ Tiny Tim/Spitzer, developed by John Kiist for the Spitzer Science Center. 
The Center is managed by the California Institute of Technology under a 
contract with NASA 

^ SPIRE Operating Manual: http://herschel.esac.e3a.int/ 
hcss-doc-5 . 0/print/spire_um/spire_um. pdf 



Figure 4. Optical HST observations of Beta Leo after PSF subtraction. 
North is up, east is left. The green circle has diameter 8" or 88 AU. The 
residual halo of light discussed in the text is contained interior to the green 
circle. Saturation columns are evident to the upper right and lower left of the 
star The bright circular structure above the saturation columns is a different 
occulting spot in the HRC focal plane. 




Figure 5. SPIRE observations of /3 Leo at 250 |J.m. The disc is detected but 
not resolved at 25()|xm. At 350 there is a la detection, giving a 3 cr upper 
limit on the disc and star emission of 2 1 mJy and the non-detection at 500 
|a.m also gives an 3 ct upper liinit of 12.3 mJy, which is consistent with the 
fluxes expected from the SED modelling. 




RA Offset, arcsec 

Intensity 



22.5"±1.2x 18.7"±1.1 with a PA of 118°, consistent with the PA 
of the resolved disc (125±15°) from the PACS images. Compari- 
son of the disc size with that of the SPIRE PSF (Neptune), found 
no significant extension at 250 \im. 

2.8 Interferometry 

Akeson et al. (2009) present interferometric observations of /? Leo 
at 2 [xm using the FLUOR instrument at the CHARA interferom- 
eter. These data show a short-baseline visibility deficit, indicating 
that the source is somewhat extended. Akeson et al. (2009) sug- 
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gest a model in which the inferred few percent MR (Near Infra- 
red) excess is due to a population of dust grains within the field 
of view (<4.6 AU), comprised of a population of blackbody grains 
that extends inward to the sublimation radius (0.12 AU for /3 Leo 
for a sublimation temperature of 1600K), which maximises the 2 
[im emission. This places a lower limit on the fractional dust lumi- 
nosity (defined as L/ij/L* i.e. a measure of the IR excess for each 
component.) for this hot population of 8.0±l.lx 10~^, which is 
larger than the value of 2.7 x 10^^ for the IRS emission (Chen et al. 
2006). The strongest spatial constraint from the interferometry is 
that there is dust with a fractional luminosity of 8 x 10^'' within 4.6 
AU. 

Stock et al. (2010) present new 10 |J.m nulling interferometry 
data from the Keck Interferometer NuUer (KIN), which suggest that 
there is no significant resolved emission at 10 |J.m in the KIN beam 
(FWHM of 0.50" X 0.44"(Colavita et al. 2009)). Limited by its 
small field of view (0.6"), these observations are unable to detect 
extended structures beyond a radius of 3 AU at the distance of /? 
Leo (11 pc). Moreover, due to its complex transmission function, 
KIN observations are most sensitive to extended structures with 
radii from 0.1-1 AU. Stock et al. (2010) also present an N-band 
(8-13 \va\ ) nulling interferometry detection with BLINC of a null 
of 1.74±0.3%, which leads to a model-based estimate of an excess 
flux of 250-I-/-50 mJy. BLINC is sensitive to the region between 1 
and 9 AU (between 0.12 and 0.8"). Stock et al. (2010) quote limits 
for various uniform disc models: 110 mJy for a disc between and 
I AU, 240 mJy for a disc from 0-2 AU and 370 mJy for a 1 to 2 
AU ring. Stock et al (2010)'s preferred model is a ring from 2-3 AU 
with a flat surface density profile, and a flux of 250±50 mJy at 10 
(im. This flux is used as a constraint in the SED modelling. 



and were then combined linearly with equal weight to come to a 
final chi-squared (Xr„, = Xwineias +Xw»ne35)- 

The emission spectrum of the photosphere of /? Leo was dis- 
cussed in § 2.1. To calculate the emission from the dust grains 
in the model they were assumed to have a size distribution with 
n{D) (X D^'^'^\ where D is the grain diameter, which is the stan- 
dard solution for a theoretical collisional cascade (see Dohnanyi 
1969), that is truncated at a minimum and maximum grain size. 
The minimum grain size was treated as a free parameter with a 
range between 0.01-100 times Dm, where Dbi is the largest grain 
size blown out of the system by radiation pressure and depends on 
grain composition and stellar properties. The maximum size was 
fixed at 1cm since larger grains have a negligible contribution to 
flux from a model with this size distribution. The grains were as- 
sumed to have a silicate core (amorphous olivine) and an accreted 
mantle of organic refractories produced by UV photoprocessing of 
ice (as used in Li & Greenberg (1997); Augereau et al. (1999)). A 
range of compositions was also tried, with amorphous silicate frac- 
tions varying from 0% to 90% by volume and with porosities (i.e. 
vacuum fraction of grain by volume) from % to 95 %. Dielec- 
tric constants were calculated from tabulated laboratory values (Li 
& Greenberg 1997; Augereau et al. 1999) using Maxwell-Garnett 
effective medium theory. The optical properties of the grains were 
calculated using Mie theory, Rayleigh-Gans theory and Geometric 
Optics in the appropriate size regimes (Bohren & Huffman 1983). 
This composition was the best fit found by minimising x^ across 
a grid representing the possible compositions. However, we do not 
give constraints on the composition because the model used con- 
tains uncertainties in the calculation of the optical properties and 
we do not want to overempahise this aspect of the modelling. 



3 MODELLING 

To determine constraints on the radial distribution of the emission 
seen in the resolved 100 and 160 ]im images we first considered a 
model of the disc structure. For example in §3.1 we consider the 
simplest possible model. This model is composed of a single ax- 
isymmetric disc component, defined by four free parameters: inner 
radius (ii„), outer radius (rout), inclination (i) and surface density 
distribution E, which is assumed to have the form S ex r^ . The disc 
opening angle (which sets the disc height) is assumed to be 5° but 
this parameter is unconstrained by the modelling process. The flux 
from annuli in the disc at different radii was determined assuming a 
grain composition and size distribution that were constrained using 
the emission spectrum. 

In order to constrain the dust location, the width of the disc, 
the inclination and the surface density profile, a grid of models was 
run. The data used in this modelling is listed in Table 2, which in- 
dicates which observations were used to constrain the SED and the 
images. The model images were convolved with a PSF (image of a 
Boo taken in the same observing mode as the observations of (3 Leo 
and reduced in the same way) and compared with the observations 
using the images as well as the linecuts both in the direction of ex- 
tension (PA 125°) and perpendicular to the extension (PA 35°). The 
use of both linecuts allowed us to constrain the inclination simulta- 
neously with the radial morphology. Both wavelengths were fit si- 
multaneously, and a joint best fit was determined by minimising the 
combined reduced x^. The reduced x^ of the fit (Xr = ^~^° ) 
where v is the number of free parameters and 1 represents a good 
fit) to each of these six pieces of observational data were calculated 



These models were fit to the IRS spectrum, BLINC point, 
MIPS points, PACS points and SCUBA-2 upper hmits. Table 2 lists 
the data used in the modelling and indicates if it was used to con- 
strain the SED, images or both. As the IRS spectrum has no obvi- 
ous features to fit to, to calculate Xsed we chose 6 windows in the 
spectrum, each ~l-2 \an wide, with a constant signal to noise ratio 
in the window. This gives 11 fluxes (6 IRS, 1 BLINC, 2 MIPS, 2 
PACS) and 2 upper limits with which to calculate the Xsed- 

The IRS spectrum is ~2% lower than the MIPS observation at 
24 [xm as the MIPS observation is marginally resolved, as discussed 
in § 2.5. When calculating the MIPS 24 |j.m photometry a larger 
aperture was used to ensure all the flux was included, whereas the 
IRS spectrum assumes a point source. A synthetic model image 
at 24 (j.m was compared with the MIPS image to check that the 
model reproduces the marginal extension observed. The total flux 
was compared to both the IRS and MIPS points and the ~2% dis- 
crepancy is a small fraction of the total flux and does not affect 
any conclusion about the goodness of fit of the models. There is a 
similar discrepancy between the BLINC flux at 10 (im and the IRS 
spectrum, which is again small enough not to affect the conclusions 
of the modelling within the uncertainties. 

We considered models of increasing complexity. The initial 
model (described in § 3.1) is a 2 component model with a hot inner 
disc and a cold outer disc that has been resolved with PACS. The 
second model has 3 components, a hot inner disc, a warm com- 
ponent and a cold outer disc and is described in § 3.2. The third 
model consists of a single eccentric planetesimal population and is 
described in 5 3.3. 
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Table 2. Table listing the data used in the modelling process, indicating if it was used to constrain the SED or the images or both. Where the table indicates 
that the observations were used as a surface brightness constraint this indicates that model images were produced at this wavelength to check that the predicted 
model surface brightness was compatible with the observed limits. 



WavelengthAVaveband 


Flux (Jy) 


Error 


Instiiunent/Survey 


Resolved 


Constraint on Modelling 


Reference 


1.6 \i.\ji 


N/A 


N/A 


HST ACS 


No 


Image Surface Brightness (Upper Limit) 


§2.6 


3.6 \im 


49.435 


1.085 


Spitzer IRAC 


No 


SED 


Stock etal. (2010) 


4.5 i^m 


30.971 


0.693 


Spitzer IRAC 


No 


SED 


Stock etal. (2010) 


5.8 \xm 


20.450 


0.458 


Spitzer IRAC 


No 


SED 


Stock etal. (2010) 


6.75 |J.m 


15.07 


0.081 


Spitzer IRS 


No 


SED 


Chen et al. (2006) 


8.0 [xm 


10.759 


0.244 


Spitzer IRAC 


No 


SED 


Stock etal. (2010) 


9.0 


8.094 


0.096 


AKARI 


No 


SED 


Murakami et al. (2007) 


9.5 i^m 


8.22 


0.082 


Spitzer IRS 


No 


SED 


Chen et al. (2006) 


10.5 ^m 


6.55 


0.05 


BLINC 


Yes 


SED 


Stock etal. (2010) 


11.3nm(Nband) 


5.82 


0.48 


Gemini MICHELLE 


No 


Surface Brightness limits, SED 


§2.3 


13. nm 


4.29 


0.024 


Spitzer IRS 


No 


SED 


Chen et al. (2006) 


15.5 nm 


3.36 


0.073 


Spitzer IRS 


No 


SED 


Chen et al. (2006) 


18.1 iJ.m(QBand) 


2.36 


0.31 


Gemini MICHELLE 


No 


Surface Brightness limits, SED 


§2.3 


20.0 ^m 


2.30 


0.021 


Spitzer IRS 


No 


SED 


Chen et al. (2006) 


lAyim 


1.647 


0.033 


Spitzer MIPS 


Marginally 


SED, Image 


Stock etal. (2010) 


29. nm 


1.68 


0.009 


Spitzer IRS 


No 


SED 


Chen et al. (2006) 


33.5 ^m 


1.54 


0.014 


Spitzer IRS 


No 


SED 


Chen et al. (2006) 


lOyim 


0.743 


0.052 


Spitzer MIPS 


No 


SED 


Stock etal. (2010) 


100|j.m 


0.480 


0.030 


Herschel PACS 


Yes 


SED, Image 


§2.2 


160 urn 


0.215 


0.032 


Herschel PACS 


Yes 


SED, Image 


§2.2 


250 i^m 


0.051 


0.012 


Herschel SPIRE 


No 


SED, Image 


§2.7 


350 i^m 


<0.039 


N/A 


Herschel SPIRE 


No 


SED upper limit 


§2.7 


450|j.m 


<0.050 


N/A 


SCUBA-2 


No 


SED upper limit 


§2.4 


500|j.m 


<0.015 


N/A 


Herschel SPIRE 


No 


SED upper limit 


§2.7 


850 i^m 


<0.006 


N/A 


SCUBA-2 


No 


SED upper limit 


§2.4 


870 i^m 


<0.20 


N/A 


HHT Sub-mm Observatory 


No 


SED upper limit 


Holmes et al. (2003) 



Figure 6. Spectral Energy Distribution (SED) of /3 Leo. The photosphere of (3 Leo is fitted with a Kurucz model profile (Lstar=14.OL0, Tstar=8660K) 
fitted to the 2MASS fluxes and shown with a dark solid line. The IRS spectrum (solid line from 5-30 |j.m) of Chen et al. 2006, MIPS fluxes (crosses), BLINC 
lO.l^m (upright cross) and Gemini MICHELLE fluxes in 1" radius aperture form imaging described in§2.2 (diamonds) and PACS fluxes in 20" aperture from 
imaging described in §1 (squares). Upper limits are 3- u from SPIRE (350, 500 i^m), SCUBA-2 (450 and 850 ^m), SCUBA (higher 450 |xm, 850 |xm values) 
and Bolocam (1.1mm), The excess is fitted with a 2 component realistic grain fit described in §3.2 (dark solid line - hot component is shown with a dashed 
line and cold component is shown with a dotted line). 
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Figure 7. Model PACS images for best fit 2 component model - a cold disc from 15-70 AU with a surface density profile S oc r~^'^ and an unresolved hot 
component at 2 AU. This is shown to the same colour scale and pixel size {l"at 100 and 2" at 160 (xm) as the observations shown in Figure 1. The 100 |j.m 
image is on the left, the 160 |J.m is on the right. The residuals (observations -model) are below, with 100 on the left and 160 ^.m on the right respectively. 




RA Offset, arcsec 

Intensity 



D.0S40 ([inear) 3.ago4>3 



RA Offset, arcsec 

Intensity 




Offset, arcsecs 
Intensity Scale 



5.98eJ (lineai-l 9.S7eJ 




3.1 2 Component Model 

The ranges of model parameters tested for the 2 component model 
were i?i„:5 to 80 AU (1 AU intervals), Rout- 20 to 150 AU (5 AU 
intervals), inclination: 0° to 90° (5°) intervals, where 0° is edge 
on), surface density index 7: to -3.0 (0.5 intervals). The values 
for gamma were chosen to cover possibilities such as the surface 
density distribution expected from grains being blown out of the 
system by radiation pressure (7=- 1.0), and that of the Minimum 
Mass Solar Nebula (MMSN) (7=- 1.5). The best fit 2 component 
model to the PACS 100 and 160 (xm images, when considered iter- 
atively with the SED fitting described in §3 was found to be a ring 
between 15±10 to 70±5 AU, with an inclination of 55±5° from 
edge on, and a surface density profile index 7=-1.5±0.5. minimum 
grain diameter of 0.5 x Dm (3 [im), a fixed maximum grain size of 
1cm and a composition with a silicate fraction of 20 %, a porosity 
of 20 % with the rest of the grain composed of organic refractories, 
with no ices present. The hot dust is assumed to lie from 2-3 AU 
with a flat surface density profile with a silicate fraction of 60% and 
a porosity of 20% and a minimum grain diameter of 0.6 fi.m. 

The fractional luminosity of the cold component is 3x 10~^. 



As there are no obvious features in the IRS spectrum, the main 
constraint on the composition and size distribution comes from al- 
lowing the appropriate range of temperatures to be present, given 
the constrains that the grains are in the region 15-70 AU; blackbody 
grains would have to be at a radius of 19 AU to achieve the observed 
temperature of 120 K. Although our model provides a consistent fit 
to both image and SED, it is not expected that the composition has 
been uniquely constrained by this process. 



To evaluate the fit of the model to both the SED and the 

+ XSED, where Xrres = 

, is the combined reduced 



images, we use X?,o„b,„ed=X?res 
2 -1-2 H ^ 



duced x^ would be 1. This best fit model had Xr^ 

.2 _i Af, ^2 _i CO „„j ,,2 



=4.5 with 



Xwineioo = l-46, Xr!inei6o = l-58 and X5BD=l-46. For a perfect fit, 
x'i u- J would be 3. The best fit y^ ^. ^ model image is 
shown on the left of Figure 7 and can be compared directly with 
the observed disc (Figure 1) which has the same colour scale. 

There was a difference between the models preferred by the 
SED and the images. The SED best fit requires a closer inner 
edge (ri„=10 AU) to produce enough flux in the 10-70 iJ.m region 
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whereas the images prefer a shghtly larger inner radius (ri„=25 
AU) with the same size distribution and grain composition.The best 
fit value of Tin is therefore a compromise between these two val- 
ues. To constrain the error on the inner radius, a Bayesian infer- 
ence method is adopted (e.g., Lay et al. 1997, Akeson et al. 2002, 
Pinte et al. 2008), in which each model is assigned a probability 
that the data are drawn from the model parameters. In cases where 
the Bayesian prior has a uniform probability distribution, as is the 
case here, this probability is T'=Paexp~^ ' ■^, where x^ is the unre- 
duced chi squared. The normalisation constant, Pq, is chosen so 
that the sum of the probabilities over all models in the grid is unity. 
Once this is done for all models in the grid, the probability distri- 
bution for a given parameter can be derived by marginalising the 8- 
dimensional probability hypercube against the other 7 dimensions. 
An example of this process for the disc inner radius is shown in 
Figure 8, which shows that the distribution for possible locations 
of the inner edge peaks at 15 AU with a range of 10-40 AU. 

It is also possible to marginalise 2 parameters against the other 
six, creating a 3D plot with the probability distribution which al- 
lows us to examine the dependence of parameters upon each other. 
Figure 8 also shows such a plot of the marginalisation for the outer 
radius and gamma, the exponent of the surface density power law. 
The outer radius and surface distribution are not well constrained 
by fitting the SED alone, partly due to the lack of constraints at 
wavelengths longer than 160 [im. The images suggest that the outer 
radius required is >60 AU, but these parameters are degenerate as 
the steeper the surface density profile, the larger the outer radius 
needed to provide a good fit. The best fit choice of rout =70 AU 
and surface density index 7=- 1.5 represent the best compromise 
between these 2 factors. The grain size constraints on composition 
and minimum grain size come solely from the SED. 

Compared to the two-component disc model presented in 
Stock et al. (2010), which has a ring of 1 \x.m diameter carbona- 
ceous grains in a ring at 2-3 AU and a disc from 5-55 AU with a 
flat surface density profile (7=0) composed of silicate grains with 
radii from 5-1000 \an. The Stock et al. (2010) model is driven by 
the resolved structure by BLINC and the excess level of 250 +1- 
50 mjy. Constraining from the excess range of 5.8-10.5 \iva region 
(IRAC 5.8 and 8 |J.m, BLINC 10.5 ^m), a temperature of 600 K is 
imposed for the inner hot component, resulting in an excess peak- 
ing at 5 11m. In comparison, the hot component in our model has 
a lower temperature due to a different size distribution and grain 
composition which means that it matches the Stock et al. (2010) 
model well at >5 jim but was not fitted in the IRAC bands due to 
differences in the photospheric model. There is a difference of 5% 
in the photospheric fit used in the IRAC bands, but the photospheric 
subtraction at 10 and 24 \im produces values that are within the er- 
rors of those used by Stock et al. (2010), despite the difference in 
the stellar fit. When constraining the cold and warm components, 
we primarily consider the excess as >5 \an and refer to the Stock et 
al. (2010) model for a more thorough treatment of the hot excess. 

The total mass in the collisional cascade, with the assumed 
size distribution of n(D) oc D^'^'', scales with Mtot oc \/Dmax- 
Thus scaling to the maximum grain size of 1 cm gives a dust mass 
of 2.1 X 10"" M®, with 3.2x 10"* M® of this in the hot compo- 
nent and the rest in the cold component. Assuming the collisional 
cascade extends up to bodies of Dmax=^ km then the total mass 
would be 1.2 M®. PSF variation can have an important effect on 
the observed residuals, so we repeated the modelling process using 
the Vesta PSFs and found no significant differences in the best fit 
model within the quoted errors, indicating that PSF variation has a 
negligible impact on the results. 



Figure 8. The top figure shows Bayesian Marginalisation of the best-fit 
model grid showing the probability distribution for the inner radius fit The 
bottom figure shows Bayesian Marginalisation of the best-fit model grid 
showing the probability distribution between the outer radius and the power 
law index of the surface density distribution, gamma. The highest proba- 
bility is represented by the lightest colour, with a maximum probability of 
0.53. . 
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3.1.1 Comparison with MCFOST Results 

We also compared the results of out model to the best fit found 
by the MCFOST code (Pinte et al. 2006), a Monte Carlo radia- 
tion transfer code in which a star radiates isotropically in space 
and illuminates an azimuthally symmetric parametrised disc. This 
was used to fit both the Herschel 100 and 160 (J-m images and the 
entire SED of /3 Leo. A grid of approximately a million models 
was run with the range of parameters described above and the joint 
X^ of the images and SED was minimised. The best fit model de- 
rived from fitting the SED and PACS images with MCFOST gives 
qualitatively similar best fit parameters as those derived from our 
IDL modelling suite (see § 3.2). As the two approaches give sim- 
ilar results, this validates the models within the constraints of the 
assumptions made about spatial and size distributions and compo- 
sitions. The SED was weighted more heavily in calculating the x^ 
in the MCFOST models, and the tension between the best-fit for 
the SED and images results in a compromise with the model pre- 
ferring either most of the mass in a relatively narrow ring around 
30-40 AU (as indicated by the SED), or a more extended disc start- 
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ing further in (with smaller grains) and extending further out with 
a shallower surface density profile, which provides a better fit to 
the PACS images. As our model grid weighted the goodness of fit 
to the image surface brightness profiles more heavily than the SED 
fit, we consider the more extended disc as our best fit model. The 

24 |i.m image was not included in the MCFOST grid, and the exten- 
sion seen in this image provides evidence against the narrow ring 
interpretation. 

3.2 Three Component Model 

In § 3.1 we describe a 2 component model of the j3 Leo debris disc 
system. However, there is some compromise between the best fit 
models indicated by the SED and PACS images which leads to an 
uncertainty in the location of the inner edge of the disc between 
15-30 AU, as this region is within the PACS beam size. This ambi- 
guity is due to limitations from the initial assumption of a disc with 
an inner radius, an outer radius and a continuous surface density 
profile. By considering the possibility of a more complex 3 compo- 
nent model we can place better constraints on the inner region of 
the disc. By combining the SED and surface brightness limits from 
the mid-IR and scattered light images we can place limits on pos- 
sible 3 component models for the system and constrain any warm 
dust population between the hot and cold components. This 3 com- 
ponent model has 7 parameters: radius and cross-sectional area of 
hot dust ( rfiot, Ghot), radius and cross-sectional area of warm dust 
(rmarm, o-ujarm), radlus and cross-sectional area of cold dust {Vcoid, 
Ccoid) and composition of the warm component. We fix the compo- 
sition of the hot and cold components to that inferred in §3.2, since 
we know that this results in an appropriate range of dust tempera- 
tures at a given distance to fit the PACS images. The parameters of 
hot dust are set by the limits from the lOpim interferometry - here 
we restrict the parameter space to That of 2, 5 or 8 AU. The param- 
eters of the cold dust are set by the PACS imaging as the best fit 
cold disc described in § 3.2. 

The main contribution from the cold component to the emis- 
sion at 10-70 (i.m comes from near the inner edge of the distribu- 
tion. We tried inner radii of the cold component of 15 AU, 20 AU, 

25 AU and 30 AU, all of which are consistent with the PACS imag- 
ing. The warm component can then be constrained through SED 
fitting (See Table 2), primarily considering the fit to the IRS and 
MIPS photometry, and then checked for consistency with the sur- 
face brightness limits from the 11.3 and 18.1 (i.m Gemini imaging 
(§2.3) and with the 24^m MIPS image (§2.5). 

The constraints on the warm component can be examined by 
considering the percentage of the total 18 (i.m flux that is produced 
by each component in the disc. If we assume the hot component is 
at 2 AU, then for this to fit the BLINC observations at 10 |i.m, the 
warm component needs to produce 19% of the total 18 (i.mflux, (48 
mjy). The cold component with an inner edge at 15 AU produces 
8 1 % of the 1 8 ]im flux, which means that there is no need for a 
warm component. This is shown in the first panel of Figure 13. If 
the inner edge is moved to 20 AU then only 70% of the 18 \an 
emission is produced by the cold component, meaning that warm 
emission would have to account for 10% (25 mJy) of the 18 ]im 
emission. If the inner edge is at 30 AU then the cold disc only 
accounts for 51% (129 mJy) of the total 18 (xm emission, leaving 
29% (74 mJy) that must come from the warm component. 

If we repeat this analysis with a hot component at 5 AU (which 
produces 29% of the 18 [im emission) we find that the inner edge 
of the cold component cannot be at 15 AU because this produces 
too much 18 (J-m emission, so the inner edge must be >20 AU. For 



an inner edge at 25 AU, we need a flux of 23 mJy at 18 (im (9% of 
the total 18 ]im flux) from the warm component. Similarly for an 
inner edge at 30 AU we need a flux of 5 1 mJy (20% of the total 18 
(im emission). These possible configurations are shown in Figure 
13. 

Although the above arguments imply that there are many pos- 
sible hot dust radii and cold component inner edges, there are fur- 
ther constraints on the spectrum and many require unphysical as- 
sumptions about grain properties. There are two best fits to the 
3 component model. The first has dust at 2 AU, 9 AU (Temper- 
ature of 160K) and 30-70 AU (7=-2.0 for the cold component), 
whilst the second has dust at 5 AU, 12 AU and 30-70 AU (7=- 
2.0). The same composition is used as that of the best fit 2 compo- 
nent model (See § 3.1). The inteferometric constraints give slight 
preference to the first scenario with the dust at 9 AU. This gives 



us a Xr„ 



2 



=4.45 with Xr, 



lineiQQ' 



= 1.11, 



Xrli' 



= 1.31 and 



Xs_E_D=2.03, showing that the 3 component model is a better fit 
than the 2 component model to the observations. As the 3 compo- 
nent model has more free parameters than the 2 component model 
(19 and 16 free parameters respectively) it is expected to have a 
lower x^ than the 2 component model, but as the reduced x^ takes 
into account the number of free parameters, v they can be directly 
compared. The significance of this decrease in x^ is discussed in § 
3.4. The SED with this fit is shown in Figure 11. 

3.3 Eccentric Ring Model 

Two other stars that harbour hot dust at ~1 AU resolved interfero- 
metrically are HD69830 (Beichman et al. 2006; Smith et al. 2009a), 
and 7] Corvi (Smith et al. 2009b). Both these systems can be fit with 
a single continuous planetesimal population, in the form of a very 
eccentric (e>0.9) ring (Wyatt et al. 2010) with a pericentre at the 
location of the hot dust and an apocentre corresponding to the cold 
belt for rj Corvi, and an as yet unseen cold population for HD69830. 
We used this model (Wyatt et al. 2010) as a third option to try to ex- 
plain the emission seen around /? Leo. This model assumes that the 
star formed with an eccentric planetesimal population and that this 
population has been evolving due to steady state collisional erosion 
for the 45 Myr age of the system. The model has two dust com- 
ponents: cold dust coincident with the planetesimals and hot dust 
created in collisions at pericentre that is being removed by radia- 
tion pressure. The model parameters are the pericentre, apocentre, 
maximum planetesimal diameter and the size of the grains being 
removed by radiation pressure. The best fit model was found by 
minimising the reduced x^ fit to the SED and the 12, 18, 100 and 
160 (im images across a grid of model parameters. The range of 
parameters tested to find the best fit were: pericentre from 0.5-5 
AU (interval 0.5 AU), apocentre from 40-140 AU (interval 10 AU) 
and maximum planetesimal diameter between 500-3000 km (inter- 
val 500 km). The population of grains being removed by radiation 
pressure were assumed to be a single size, and values from 20-100 
% of the blow out size were tested. 

The SED for the best fit model is shown in Figure 10. This 
model has a single planetesimal population with pericentres at 2 
AU and apocentres at 65 AU. The eccentricity was therefore 0.96 
and the maximum planetesimal diameter was 2000 km . The grains 
blown out from pericentre have a size of 2 (xm. The resulting mass- 
loss rate given by this model is 0.005 M^/Myr, giving a mass of 
0.45 M0 for the parent planetesimal population. It is therefore pos- 
sible to explain the /? Leo disc using one continuous planetesimal 
population, but the possible origins of such an eccentric planetes- 
imal population are an issue (Wyatt et al. 2010). The combined 
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Figure 9. The SED fit for the 3 component model of the /3 Leo disc. The SED shows the individual contributions (dotted lines) from a hot component at 2 AU, 
a warm component at 9 AU and a cold component at 30-70 AU. The stellar photosphere is shown as a straight line and the cumulative disc distribution is a 
solid black line. 
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reduced y^ for this model is y^ ^. =4.47, showing that the 
eccentric ring model has a slightly lower y^ ^ ^ than the 2 
component model, but is a slightly worse fit than the 3 component 
model. The significance of this change in reduced y^ is assessed in 
§3.4. The fit to the surface brightness profiles of the 2 component, 
3 component and eccentric ring models at 12, 18, 100 and 160 iJ.m 
are shown in Figure 1 1 . The 24 |j.m linecut shows that the eccen- 
tric ring model appears marginally more extended when compared 
with the observations. The FWHM of the eccentric ring model at 
this wavelength is 7."0 x 6."7 compared to a FWHM of 6."88 x 
6. "61 for the /3 Leo observations. This means that the eccentric ring 
model is slightly less favoured, but cannot be ruled out. 



3.4 Ambiguity of the Models 

This addition of an extra parameters in the 3 component and ec- 
centric ring models improves the y? ^ . ^ from 4.50 to 4.45 and 
4.47 respectively, where a perfect fit to the resolved images at 100 
and 160 um and the SED would have y? ^ ^ =3.0. However, 

r /vi combined ' 

in general adding an extra parameter will improve the fit of any 
model, so is adding this extra parameter justified? Although the re- 
duced x^ compensates somewhat for this as it is divided by the 
number of free parameters, u, when considering models with large 
numbers of free parameters ( 1 9 in the case of the 3 component case) 
we need to avoid the problem of "overfitting" or choosing a more 
complex model than is warranted by the data. We assess this using 
the Bayesian Information Criterion (BIC): 



BIC = Nln{xlc 



-klnN 



(1) 



where N is the number of data points, k is the number of free pa- 
rameters and y? ^ ^ is the minimum combined reduced y*^ for 

^^' combined '^ 

the model. The BIC considers the fit of the model but penalises the 
model for extra parameters (see Wahhaj et al. 2005; Liddle 2007). 
It penalises the model more strongly for extra parameters than the 
reduced y^ and so provides a more stringent test of the validity 



of applying more complex models. The BIC value for the 2 com- 
ponent model is 51.1, for the 3 component model it is 49.7, and 
for the eccentric ring model the BIC is 51.9. A lower value of the 
BIC is preferred; a difference of 2 between BICs of different mod- 
els indicates positive evidence against the higher BIC value, and 
a difference of 6 indicates strong evidence against the higher BIC 
value. Therefore our results indicate that there is no real prefer- 
ence between the 3 models since they all reproduce the data well. 
Furthermore it is not unreasonable to assume that a realistic disc 
system may well be more complicated that the simple parameteri- 
sations used in these models. 

These three models all provide a very good fit to the observa- 
tions, but further data could allow us to resolve this ambiguity. The 
best way to distinguish between the possible 3 component models 
and the 2 component model would be to resolve the inner edge or 
the warm component. This may be possible with very deep mid-IR 
imaging in the N band. Resolving the disc at this wavelength would 
give different surface brightness profiles for the different models 
as shown in Figure 12. The Spitzer image at 24 (j.m already sug- 
gests that the extent of the eccentric ring model may be too large. 
A quadratic subtraction of the FWHM of the beam from the image 
FWHM diameter (after PSF subtraction) indicates a diameter of ~ 
4 " FWHM for the disc at this wavelength, significantly less than 
the value indicated in Figure 13. 

It may also be possible to resolve the inner edge of the cold 
disc in scattered light. We created simulated ACS images for all 
three models to compare to §2.6. Mie theory was used to calcu- 
late the albedo of grains in the model and a bandpass of 0.45-0.72 
(j.m was assumed for the F606w filter. A Henyey-Greenstein phase 
function (Henyey & Greenstein 1941) was used to approximate the 
asymmetric scattering by small particles with an asymmetry param- 
eter (g) of 0.3 corresponding to slightly forward scattering. Due to 
uncertainties in scattering properties we might expect our predicted 
surface brightnesses to have significant uncertainties. However, the 
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Figure 10. The SED of /3 Leo modelled using a single planetesimal population with pericentre at 2 AU and apocentre at 65 AU. The dashed, dotted and 
solid lines correspond to the contribution from the collisional cascade, blow out grains and total emission spectrum respectively. The diagonal solid line is the 
stellar spectrum, the Spitzer IRS spectrum from 10-40 i^m is plotted with a solid line, the asterisks are the MIPS points, the triangles photometry from Gemini 
Michelle and the squares PACS photometry. These are all shown with photospheric subtraction and have been colour corrected where necessary. 



3 




10.0 100,0 

Wavelength (jjrn) 



1 000.0 



Figure 11. Surface brightness profiles for the /3 Leo disc at 12 |j.m (top left), 18 |xm (top centre), 24 i^m (top right), 100 |xm (bottom left) and 160 \in\ (bottom 
right) showing the fit of all 3 models to the observations. These are line cuts taken parallel to the major axis of the disc (125°), summed over a width of 4 
pixels. The observed profile is shown with a solid line, the model profiles with coloured dashed lines (red for the 2 component model, blue for the 3 component 
model and green for the eccentric ring model). The bottom panel of each plot show the residuals i.e. the observed linecut minus the model linecut to show the 
goodness of fit of the model. 
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Figure 12. Predicted 1 1.3 |xm surface brightness profiles for the three different models of the (3 Leo disc: the 2 component model with a hot disc at 2-3 AU 
and a cold disc from 15-70 AU; a 3 component model with a hot disc at 2-3 AU; a warm ring at 9 AU and a cold disc from 30-70 AU and the eccentric ring 
model. 
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Figure 13. An illustration of the possible structures of the /3 Leo disc based on the 3 different possible models of the disc. The top panel shows the 2 component 
model, which has belt at 2 AU, a gap from 2-15 AU and an extended disc from 15-70 AU. The middle panel shows the 3 component model which has a belts 
at 2 AU, 9 AU and 30-70 AU. The bottom panel shows the eccentric ring model as described in the text. 
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morphology for the emission (and relative level of surface bright- 
ness for different models) should be correct. 

The ACS coronograph has an inner working angle of 1.8 " so 
in all the models the hot component is hidden behind the occult- 
ing spot. For the 2 component model the peak surface brightness in 
scattered light would be a radius of 1.35 " (i.e. behind the occulting 
spot) with a surface brightness of 1.7 x 10^'^ mJy/pixel, assuming 
a pixel size of 25 mas. Just outside the inner working angle of the 
coronograph (1.9 " radius) the disc would have a surface brightness 
of 1.3x 10~^ mJy/pixel and fall off oc r"''"''. The three component 
model has peak surface brightness of 2.1 x 10"'* mJy/pixel at a ra- 
dius of 2.7 " and the eccentric ring model has a peak surface bright- 
ness of 4.6x 10^" at a radius of 3.6". In an annulus 2.8" - 3.0" 
the residual halo in the HST ACS observations shown in §2.6 has a 
median brightness 5.3 x 10"'* mJy/pixel, which is consistent with 
the three component model. The halo brightness decreases steeply 



to the sensitivity-limited sky value at 4.0" radius. If the halo were 
due to grain scattering, we would expect a radial brightness profile 
proportional to r^^ if the star was embedded in a uniform density 
sheet of dust (7=0). The measured radial brightness profile is pro- 
portional to r~^, implying that 7~-3. This is comparable to the 
dust radial profile for the southwest midplane of the /3 Pictoris dust 
disc (Kalas & Jewitt 1995; Golimowski et al. 2006). However, this 
is significantly steeper than 7=-2.0 for the three component model. 

The nature of this residual emission is unclear - it is difficult 
to find a disc model that can both produce sufficient scattering in 
this region without producing a very bad fit to the SED. Models 
with an outer edge at 44 AU, as seen in the ACS observations or 
with 7<-2.0 to mimic the steep drop-off seen are also a very bad fit 
to the extended emission seen with PACS at 100 (im and 160 |J.m. 
However, it could be reconciled with the two component model if 
we are seeing a separate small grain population in the optical ob- 
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servations that has a steeper surface density profile. If this emission 
was real, then it could be used to place limits on the albedo of the 
grains in the eccentric ring model, or possibly rule out this scenario. 
Currently, the observations favour the two component model as no 
peak is seen at 2.7 or 3.6". Further optical imaging could confirm 
the nature and location of this emission and possibly rule out the 
eccentric ring model, which has the largest surface brightness in 
this region. 



4 DISCUSSION 

The /3 Leo debris disc system consists of multiple dust populations. 
We can find a best fit solution with either 1, 2 or 3 components. 
These are shown in Figure 13. The 2 component model is the sim- 
plest and has dust at 2-3 AU and 15-70 AU. There are two best 
fits to the 3 component model. The first has dust at 2 AU, 9 AU 
and 30-70 AU, whilst the second has dust at 5 AU, 12 AU and 30- 
70 AU. The inteferometric constraints give slight preference to the 
first scenario. The key constraint on these models is that the warm 
(middle) component cannot be located inside an 8 AU radius with- 
out violating the 10 \im interferometry constraints. A third possible 
scenario for this system is that the emission results from a single 
very eccentric (0.96) planetesimal population. This model consists 
of a steady-state coUisional cascade with an apocentre coinciding 
with the location of the cold emission seen in the PACS images, 
with an additional population of hot, small grains created at col- 
lisions near pericentre which are subsequently removed from the 
system due to radiation pressure. 

These inferred structures give us clues to the underlying dy- 
namics of the system, as the belts are confined and separated by 
gaps suggesting dynamical interactions. The 2 component model 
has a gap between 2 and 15 AU which could contain planets, as 
could the gaps between 2 and 9 AU and between 10 and 30 AU in 
the 3 component model. 



4.1 Origin of the Hot Dust 

The origin of the hot, compact dust populations is still unclear. 
There are five possible scenarios: i) the small dust grains produced 
by collisions in the cold belt could drift towards the inner region 
due to PR drag; ii) the hot dust could be evidence of a planetes- 
imal population in a steady state collisional cascade at a few AU 
from the star; iii) the dust could be transient, produced in a mas- 
sive collision between planetesimals; iv) the system could be un- 
dergoing a dynamical instability resulting in the planetesimals be- 
ing thrown in from the outer belt and producing hot dust as in the 
Late Heavy Bombardment (Gomes et al. 2005); or v) the grains 
may be produced through evaporation of comets originating from 
the cold planetesimal belt imaged in the sub-mm, like the Zodiacal 
Cloud (Nesvomy et al. 2010). 

The first scenario concerning PR drag as an origin for the hot 
dust is improbable due to the long timescale for PR drag com- 
pared to the shorter collisional timescale for the outer disc (Wy- 
att et al. 2007). This mechanism is unlikely to produce the amount 
of dust observed in /? Leo's inner regions. Around late type stars 
radiation pressure is inefficient at removing small grains, so these 
can be transported into the inner regions of the disc via Poynting- 
Robertson drag and stellar wind drag, even in discs with optical 
depths considered too high for discs to be transport-dominated (r 
^10^* Reidemeister et al. 2010). This is because close to the star 
larger grains are preferentially coUisionally eliminated, causing a 



break in the size distribution at a critical size, below which grains 
are transport dominated. Since j3 Leo is an A type star, however, the 
larger radiation pressure blow-out size will mean that all the grains 
remaining in the disc are larger than the critical size and hence are 
coUisionally-dominated. 

The feasibility of a steady-state belt as an explanation for the 
hot dust can be examined by considering the maximum steady state 
fractional luminosity of the belt. The fractional luminosity of the 
hot component from our model is 8.1x10^'. Using equation 18 
from Wyatt (2008) for the maximum fractional luminosity of a belt 
at a given age; 



/ma. = 0.58 X 10-V"'''(dr/r)D°-'Q^'^'"'e 



-5/3 A/r-5/6,- -0.5,-1 
''age 



and assuming a radius (r) for the hot component of 2.5 AU, a width 
(dr) of 1 AU, a maximum planetesimal diameter {Dc) of 60 km, 
planetesimal strength (QJ) of 150 J/kg (the fiducial value from Wy- 
att (2008)), an eccentricity (e) of 0.05, stellar mass (M*) of 2.1 M©, 
a luminosity of 14.0 L© and an age of 45 Myr, gives a maximum 
fractional luminosity of 9.9x 10~^. For the hot dust population to 
be considered transient, Wyatt et al. (2007) concluded that the frac- 
tional luminosity must be 1000 Fmax, so /? Leo's hot dust is con- 
sistent with a steady state model, but only marginally so, and to be 
a steady state phenomenon would have to have properties that are 
significantly different from those for other A stars (e.g. unusually 
strong or large planetesimals). This implies that the third scenario 
is unlikely as again, a very large belt mass would be required due 
to the low probability of a massive collision having occurred and 
the short lifetime of the dust produced. 

The fourth scenario also interprets the hot dust as a transient 
phenomenon by suggesting that the system is undergoing major 
dynamical perturbations in a scenario analogous to the Late Heavy 
Bombardment (LHB) in the solar system. The LHB occurred 700 
Myr after the planets formed and led to a spike in the production 
of dust as comets were thrown in towards the star by a dynami- 
cal instability caused by the migration of the giant planets (Gomes 
et al. 2005). It is not known how common such events are, although 
Booth et al. (2009) estimate that they occur around <12% of sun- 
like stars. Since an LHB-like event could occur at any point in a 
disc's lifetime and as only a small percentage of stars undergo such 
an event, it is improbable that we are observing such a disturbance 
in the /? Leo disc. 

This leaves the final scenario, that the grains may be produced 
through disintegration of comets originating from the cold plan- 
etesimal belt imaged in the sub-mm. Since this requires neither a 
very high initial mass for the disc nor a low-probability event such 
as a collision or a LHB like instability, it seems more likely al- 
though the dynamics of scattering and physics of dust production 
remain poorly understood. If the hot dust does originate from a 
population of scattered planetesimals then the amount of hot dust 
could be indicative of the configuration of the planetary system re- 
sponsible for the scattering e.g. the effect of Jupiter on scattering 
of planetesimals in the Solar System (Horner & Jones 2008, 2009; 
Homer et al. 2010), This implies that there could be multiple plan- 
ets between the dust belts responsible for scattering planetesimals 
into the inner regions of the system, and so replenishing the hot 
dust population. 



4.2 Placing /3 Leo into Context 

There are currently 16 debris discs around A type stars that have 
at least one component resolved in thermal emission. These are 
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Table 3. Table comparing radius, age and fractional luminosity of debris discs resolved in thermal emission around A stars 
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160(411 








111] 






160(41] 


grains out to 
>1000 AU(4i] 


/3 UMa 


11.2(421, 


AlV 


50 


50 


1.4x10-5 


1.09±0.01 


10-5(42] 


50(301, 




(HD95418) 


18.1(421, 
100(231 










[42] 




520(44] 




HD139006 


1121421 


AOV 


461431 


46(43] 


1.2x10-5 


2.31 ±0.04 


1x10-5(421 


314(30], 


Eclipsing Binary 




18.1(421 








[43] 


[42] 




350(44] 


(companion G5V, 
O.9M0 at 0.2 
AU) 


HD181869 


10.4(421, 
18.11421 


B8V 


N/A 


N/A 


N/A 


4.12±0.11 

[42] 


IX 10-51421 


100(44] 




HD3003 


N/A 


AOV (43] 


141411 


24(41] 


2.10x10-4 

[41] 


4-6.5(43] 


2.01x10-4(43 


50(30( 
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Name 


Wavelengtl 
Resolved 


Spectral Type 


Inner 
Radius 

(AU) 


Outer 
Radius 

(AU) 


Fractional 
Luminos- 
ity (cold 
component) 


Hot Compo- 
nent 


fractional lu- 
minosity (hot 
component) 


Age 
(Myr) 


Notes 


HD23281 


N/A 


A6Vl«l 


N/A 


N/A 


N/A 


5.4-9.51431 


3.82x10-5 

[43] 


6261301 




A Gem 
(HD56537) 


N/A 


A3Vl«l 


N/A 


N/A 


N/A 


2.1-6.11431 


5.64x10-5(43 


5601431 




HD71155 


N/A 


A()Vl«l 


44I11I 


90143) 


2.5x10-51"! 


2.2-8.2143) 


8.95x10-5 

[43] 


1691301 




HD80950 


N/A 


AOVl«l 


13.6l«l 


24.0l«) 


9.62X10--"'' 

[43] 


N/A 


N/A 


80 145] 




HD141795 


N/A 


A2ml«l 


N/A 


N/A 


N/A 


4.6-6.11431 


4.43x10-5143 


4501301 





References: l)Kalas et al. (2005), 2)Holland et al. (1998) 3)Holland et al. (2003) 4)Absil et al. (2009) 5)Di Folco et al. (2004) 6)Kalas et al. (2008) 7)Su 
et al. (2005) 8)Sibthorpe et al. (2010) 9)Wilner et al. (2002) 10)Absil et al. (2006) 1 l)Wyatt et al. (2007) 12)Schneider et al. (1999) 13)Telesco et al. (2000) 

14)Moerchen et al. (2010) 15) Wyatt et al. 1999 16)Smith et al. (2009) 17)Lepine & Simon (2009) 18)Moerchen et al. (2007) 19)Kalas et al. (2000) 

20)Lagage & Pantin (1994) 21)Mouillet et al. (1997) 22)Lagrange et al. (2010) 23)Matthews et al. (2010) 24)Phillips et al. (2010) 25) This Work 26) Stock et 

al. (2010) 27)Chen et al. (2006) 28)Lachaume et al. (1999) 29)Su et al. (2008) 30)Song et al. (2001) 31)Clampin et al. (2003) 32)Weinberger et al. (1999) 

33)Silverstone et al. (1998) 35)Wyatt (2005) 36)Wahhaj et al. (2007) 37)Kalas (2005) 38)Schneider et al. (2005) 39)Moerchen et al. (2007) 40)Debes et al. 

(2009) 41)Su et al. (2009) 42)Moerchen et al. (2010)43)Smith & Wyatt (2010) 44)Rieke et al. (2005) 45)Wahhaj et al. (2003) 



listed in Table 3. Of these 16 discs 11 show evidence of multiple 
dust components, where at least one component has been imaged. 
The last 5 objects in Table 3 are A stars with evidence of multiple 
dust components including hot dust, but whose discs have not been 
resolved. Of these A stars, most are thought to be younger than /3 
Leo, with ages in the range of 10s of Myr. 

The discs summarized in Table 3 show a startling diversity, 
with radii from a few AU, i.e., C, Lep (3-8 AU Moerchen et al. 
2007) to hundreds of AU, i.e., Vega (Su et al. 2005 ). Some systems 
are confined rings with clear inner and outer edges, i.e., HR4796 
(Schneider et al. 1999) and some have very extended discs like /? 
Pic (Kalas et al. 2000). 

Figure 14 shows a plot of the fractional luminosity vs radius 
of both the hot and cold disc component for all the A stars listed 
in Table 3 as having 2 disc components. Vega, Fomalhaut and /3 
Leo are the only A stars with a hot disc component that has been 
resolved via interferometry (Absil et al. 2006, 2009; Akeson et al. 
2009). The hot dust in other discs is either inferred from SED fitting 
(as for HR4796 and HR8799) or resolved in the mid-IR (C, Lep, /3 
UMa, HD139006). This leads to potential systematic differences 
when comparing radii and fractional luminosities derived through 
different methods. These warm components could also be related 
to planetesimals near the ice-line (Morales et al. 2011). Although 
there are no clear trends in the small number of objects with known 
2 component discs, there are some interesting points to note. 

We see two main 'types' of distributions between the hot and 
cold populations. First, discs such as Vega, Fomalhaut, /3 Leo and 77 
Corvi, which have a large (>20 AU) separation between the location 
of the hot and cold dust, and whose hot dust has a high fractional 
luminosity. This leads to a steep line in Figure 14. The other type 
of disc, typified by Eta Tel, HD3003 and HR8799 show an almost 
flat or, for HR8799, declining line in Figure 14, indicating that the 
hot dust has a similar or lower fractional luminosity than the cold 
component. In the case of HR8799 the hot and cold components are 
separated by a known planetary system (Marois et al. 2008) which 
may inhibit the transfer of material from the outer to inner regions 
of the system. 

Fomalhaut, Vega and /3 Pic show the steepest increase in frac- 
tional luminosity between cold and hot components. In all 3 of 
these cases the hot component is derived from modelling fitting 
the interferometric visibility deficit seen around these objects (Ab- 



sil et al. 2006, 2009, 2010). /3 Leo shows a gradient between the 
hot and cold components, similar to that of q Corvi, an F star with 
multiple dust populations including a cold dust disc resolved with 
PACS and SCUBA at 150 AU (Matthews et al. 2010) and C, Lep, 
which has a much more compact configuration, with its 'cold' disc 
at only 10 AU. 

In the Solar System the asteroid belt (radius ~2 AU, 
f=0.8xlO"^ (Backman & Paresce 1993)) and the Kuiper Beh 
(radius~30 AU, f~10~^ (Gomes et al. 2005) would also give a 
flat line, similar to that of HR8799. The radial distribution of dust 
in these systems could be indicative of the underlying planetary 
system. HR8799 and the Solar System both have multiple planet 
systems between their dust belts which could influence the levels 
of planetesimal scattering. However, Fomalhaut and /3 Pic have at 
least one planet (Kalas et al. 2008; Lagrange et al. 2010) but still 
have high hot dust levels and a steep line on Figure 14. 



4.3 Implications for Possible Planetary Systems 

This diversity in disc structure and brightness seen in these discs 
with hot and cold components could be due to the underlying dy- 
namics of the system, which can be affected by the presence and 
distribution of planets in the system. For example, in the Solar Sys- 
tem the number of comets thrown into the inner system from the 
Kuiper belt is dependent on planetary architecture (Horner et al. 
2010). For the /? Leo disc, the radial distributions of dust that we 
have derived have implications for the architecture of possible plan- 
etary systems. For the 2 component model there is a gap in the disc 
between 2 and 15 AU, in which there could be planets that are trun- 
cating the cold disc and confining the hot inner disc. For the 3 com- 
ponent model the gaps from 3-9 AU and 9-30 AU could occur again 
due to planets truncating the discs. The size of these gaps could be 
linked to clearing by a planet due to the scattering of planetesimals 
out from the planet's chaotic zone. The area cleared increases for 
multiple planets and for planets on eccentric orbits (Bonsor et al. 
201 1). In the case of multiple planet systems these could influence 
the amount of material scattered into the inner solar system and the 
possible stable locations for dust belts, so the morphology of the 
dust belts in the Beta Leo system could indicated multiple planets 
between the cold and hot dust belts. 
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Figure 14. Plot showing the radius and fractional luminosity of the hot and cold components for the A star debiis discs listed in Table 3 as having 2 components. 
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5 CONCLUSIONS 

We have presented detailed modelling of the debris disc around the 
45 Myr old A star /3 Leo. We considered multi-wavelength data to 
construct a complete picture of this source. Resolved images taken 
at 100 and 160 iim using PACS on Herschel as part of the DEBRIS 
survey were used to place observational constraints on the radial 
location of the cold dust in this system. Resolving the cold dust is 
key to breaking the degeneracies inherent in SED modelling. We 
also use detection limits from unresolved images at 12 and 18 \im 
from MICHELLE on Gemini, at 0.6 ^m with ACS on HST, the un- 
resolved Herschel SPIRE image at 250 |i.m and detailed SED mod- 
elling including all data from the literature to gain a complete pic- 
ture of the disc. Modelling indicates that for a 2 component model 
of the system consisting of a hot and cold disc, the cold disc imaged 
with Herschel PACS lies between 15-70 AU, with a surface density 
profile E oc r"^'^ at an inclination of 55° from edge on, with the 
hot dust at 2 AU. SED fitting to observations from 5 |J.m to 1mm, 
including the IRS spectrum (Chen et al. 2006) suggest that the grain 
size distribution is consistent with that for a theoretical coUisional 
cascade with a fixed maximum grain size of 1cm, a minimum grain 
size of 3 \un (0.3 x Dm for a composition with a silicate fraction of 
20%, a porosity of 20 % with the rest of the grain composed of or- 
ganic refractories. For the hot component the minimum grain size 
is 0.6 [im (O.lxDti) and the silicate fraction was 60%. A 3 com- 
ponent model indicates that another possible configuration consist 
of an inner edge to the cold disc at 30 AU, a warm ring at 9 AU and 
hot dust at 2 AU. It is also possible to fit the observations of /3 Leo 
using a single, very eccentric (e=0.92) planetesimal population, af- 
ter Wyatt et al. (2010). This degeneracy implies that even a wealth 
of multi-wavelength data including the resolved images may not be 
enough to uniquely constrain the location of the dust when there 



are multiple populations and the edges of the belts have not been 
resolved. 

/? Leo is a similar age to Fomalhaut and Vega, which also 
have hot dust within a few AU of the star. However, its cold disc is 
smaller and closer to the hot emission that either of these systems, 
and there may be an intermediate warm component. The most anal- 
ogous disc among those resolved around A stars is Eta Tel, a 12 
Myr star with a resolved disc from 21-26 AU and hot, unresolved 
dust at <4 AU (Smith et al. 2009), but in terms of the separation 
and ratio of the hot and cold components as shown in Figure 14, /3 
Leo most resembles HD71 155. 

We have also examined the population of A stars known to 
have both hot and cold dust discs, and we see two main 'types' 
of distributions between the hot and cold populations. First, discs 
such as Vega and /? Leo which have a large (>20 AU) separation 
between the location of the hot and cold dust, and whose hot dust 
has a high fractional luminosity. The other type of disc, typified by 
Eta Tel and HR8799 show an almost flat or, for HR8799, declin- 
ing line in Figure 14, indicating that the hot dust has a similar or 
lower fractional luminosity than the cold component. In the case 
of HR8799 the hot and cold components are separated by a known 
planetary system (Marois et al. 2008) which may inhibit the trans- 
fer of material from the outer to inner regions of the system. How- 
ever, the small number of stars in this sample and the possibility of 
unknown systematic issues from comparing dust detected through 
different methods limit the conclusions that can currently be made. 
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